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Abstract. The influence of the Rashba spin-orbit coupling on the electron spin dynamics is investigated 
for a ballistic semiconductor quantum wire with a finite width. We monitor the spin evolution using the 
time-dependent Schrodinger equation. The pure spin precession characteristic of the ID limit is lost in a 
2D wire with a finite lateral width. In general, the time evolution in the latter case is characterized by 
several frequencies and a nonrigid spin motion. 

PACS. 73.21.Hb Quantum wires - 73.22.Dj Single particle states 



1 Introduction 

The control of the Rashba spin-orbit coupling is one of 
the most promising tools for the manipulation of the spin 
of the carriers within semiconductor heterostructures. In 
the last years, there have appeared several proposals of 
spin-based devices relying on this mechanism 0121 ■ Among 
them we cite the spin-FET, first proposed by Datta and 
Das P], the spin filters 0] and the spin guides [H]. 

The origin of the Rashba spin-orbit coupling in III-V 
semiconductor heterostructures lies in their asymmetry. 
The lack of space symmetry causes a local electric field 
perpendicular to the plane of the heterostructure. The 
associated relativistic correction makes the electrons feel 
an effective magnetic field which separates in energy the 
different spin states This dependence of the Rashba 
spin-orbit coupling on the interface electric field remains 
as a great advantage for spin control in semiconductors. 
Indeed, the intensity of the coupling can be controlled by 
applying a vertical electric field to the heterostructure, as 
proved experimentally by Nitta et al. [Jj. 

In this work, we shall focus on the spin precessional 
properties of conduction electrons in a quantum wire hav- 
ing a finite width. The problem of the quantum wire with 
Rashba coupling has already been studied by Moroz and 
Barnes who addressed the ballistic conductance and 
the subband structure of the wire; and by Mireles and Kir- 
czenow P|, who analyzed the spin-dependent transport 
within a tight binding model. Following a different ap- 
proach, Governale and Ziilicke [TO] have also treated the 
problem, and, more recently, Egues et al have pro- 
posed to take advantadge of the coupling of the first two 
transverse subbbands to create a spin-FET with enhanced 
spin control. In this work, we shall concentrate on how the 



coupling of the different transverse subbands, due to the 
Rashba term, affects the electron spin precession. 

The Hamiltonian representing the Rashba spin-orbit 
coupling in terms of the electron momentum (p) and the 
Paufi matrices (a's) reads 



(1) 



where represents the intensity of the spin-orbit cou- 
pling, depending on the heterostructure's vertical electric 
field. To stress the experimental feasibility we have as- 
sumed a value = 1.03 x 10~^ eV cm, as reported 
in Ref. for an Ino.53Gao.47As/Ino.52Alo.48As quantum 
well. Other experimental parameters corresponding to this 
sample are: m* = 0.05 me, e = 13.9 for InGaAs, and the 
Fermi wavevector kp — 3.5 x 10^ cm^^, corresponding 
to a 2DEG density ~ 2 x 10^^ cm^^. All calculations 
discussed below have been obtained using this parameter 
set. 

We assume a wire oriented along the y axis with con- 
finement achieved in the x direction. In the simplified 
model of a one-dimensional wire the terms containing x 
or px are dropped and the remaining Hamiltonian 
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is compatible with the observables {py:Crx}- That is, the 
energy eigenstates are characterized by their x spin com- 
ponent and y momentum (hky). In this system, the state 
of an electron having well-defined momentum and with 
arbitrary spin orientation decomposes into a combination 
of the spin eigenstates \sx+) and \sx—), leading to a spin 
precessional motion about the x axis, i.e., in the Sy-Sz 
plane. The spin evolution is characterized by an angu- 
lar frequency ujp, corresponding to the energy difference 
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between the two spin eigenstates, hujp — 2Xjiky. This pre- 
cession frequency leads to the well-known formula for the 
angular modulation of spin orientation Off = 2m* Xj^L/Tt? , 
for a channel length L. 

When a two-dimensional free motion of the electrons 
is considered, the eigenstates of the Hamiltonian 



'H2D 



2m* 



^ [PyCTx - PxCTy) (3) 



have again well-defined momentum tik and spin orienta- 
tion, but the latter depends now on the wavevector k. 
Using spinorial notation in the usual basis of |-|-) and |— ) 
states for the Sz spin component, the wave functions and 
energies for a given k read 
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(4) 
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where (j> = a,ta,n{ky / kx) is the polar angle for k. Note that 
the spinors I0J are eigenstates of a spin component orthog- 
onal to k, i.e., the spin orientation depends on the spatial 
wavevector, always pointing in perpendicular direction to 
k. As in the ID case, a state having well-defined momen- 
tum ?ik and arbitrary spin orientation decomposes into 
a linear combination of the above 2D eigenspinors Xk±- 
Therefore, the evolution under TI2D leads again to a spin 
precession of frequency ujp given by the energy difference 
between the corresponding 'up' and 'down' eigenstates. 
Namely, Tiujp = 2Ai?|k|; the same result of the ID case. 
An important difference is, however, that the plane of spin 
precession depends in 2D on the spatial orientation of k. 
Having analyzed the bulk limits, the question that arises 
is how the above pure spin precession is modified when 
the complete translational invariance of the homogeneous 
ID and 2D systems is lost in a wire with a finite width. 



2 The Model 

In order to model the quantum wire we maintain the 
translational invariance in the longitudinal coordinate (y) 
while, in the transversal direction (x) this symmetry is 
broken by a parabolic confinement 
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where luq determines the width of the wire. The Hamilto- 
nian representing the quantum wire reads 



1 * 2 2 f \ ln\ 

— + -m OJqX + — {py(Tx - PxCTy) ■ (7) 



As our system still retains one spatial symmetry we can 
reduce the spinorial 2D problem to a ID one by factorizing 
out the y dependence. 



(8) 



The continuum index ky in Eq. JSJ labels the state 
of free motion in longitudinal direction while the discrete 
index n accounts for the different transversal subbands of 
the confining dimension. Introducing this general form of 
energy eigenspinor we obtain the reduced Hamiltonian for 
the spinorial transverse dependence. 
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The transverse Hamiltonian, Eq. ((UJ, can be analyzed 
by seeking direct solutions to the stationary Schrodinger 
equation or, equivalently, solving the time-dependent Schrod- 
inger equation for given initial conditions. We shall use the 
latter approach since it yields a suitable scenario to ad- 
dress the dynamical spin evolution. For a given wavevec- 
tor ky the coupled time-dependent equations for the two 
spin components of any spinorial wavepacket: fky^{x) and 
fkyi{x) read 
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where 



Afl [ky - dx) 



Ix) 

Note that in TYq ' we have dropped the longitudinal kinetic 
contribution, not important for time evolution. 

It is well known that in quantum mechanics time evo- 
lution can provide complete information about the sta- 
tionary states of a Hamiltonian through spectral analysis. 
Indeed, for a time-independent Hamiltonian an arbitrary 
wavepacket can be expanded as a series in the stationary 
energy eigenstates 
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Using the fast Fourier transform (EFT) we can extract the 
harmonic frequencies {Snky energies) from the evolution of 
an arbitrary wavepacket. Furthermore, if this analysis is 
repeated at different spatial points the space components 
of the energy eigenspinor (pnky t (x) and (/fuky i (x) , can also 
be obtained. 



3 Results 

The Hamiltonian of Eq. Q, apart from kinetic and po- 
tential terms, contains the longitudinal and transversal 
spin-orbit contributions coupling x and y spatial motions 
with the electron spin. It is worth noticing that there 
is an analytical limit to the solution of Eq. (|10|l when 
Pxfky ^ ^kyfky. In this case the solutions are similar to 
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those discussed above for T-Cid, i-e., the eigenstates are still 
plane waves in y direction with spin oriented along x but 
now with an harmonic oscillator transverse profile. The 
corresponding energy spectrum is also analytically known 
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where s = +1 for +x spin orientation (up) and s = — 1 
for —X (down). If we substract the longitudinal kinetic 
term h A;^/2m*, common to both spin orientations, the 
remaining band structure consists of doubly split trans- 
verse subbands, linear in ky and having a common origin 
at {n + l/2)hiOQ. Another analytical limit can be obtained 
neglecting XuPyUx in iQ; in this case, spin-orbit interac- 
tion only produces a little constant shift in the energy 
levels that does not depend on spin orientation. 

As a check of our numerical procedure. Fig. 1 shows the 
results for a simulation using the above simplified Hamil- 
tonian (neglecting \RPx<yy in Eq. lO). The parameters 
used are: Uloq — 3.5 meV, ky = 1.0 x 10^ cm^^ and the 
above mentioned spin-orbit intensity. The initial input for 
the time-evolution is a Gaussian-shaped wavepacket hav- 
ing spin oriented in +z. 
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we ensure that the ini- 



Using a finite xq and a i 

tial wavepacket is composed of several transversal eigen- 
states in the low energy region. Figure 1 shows the spec- 
trum corresponding to the wavepacket 's oscillations at 
an arbitrary x point, as well as the transversal densities 
extracted for the two lowest peaks. An excellent agree- 
ment between numerical and analytical energies and eigen- 
spinors is found. Figure 2 represents the dynamical spin 
evolution for the above initial spinor, corresponding to a 
pure precession in excellent agreement with the analytical 
frequency fiojp = 2Xjiky. 

Note that, in this analytical regime, spin precession is 
independent of the transversal state, as the energy gaps 
are common to the different subbands. This statement im- 
plies that spin evolution is independent of the transver- 
sal profile of the injected particles. Numerically, we have 
checked that arbitrary initial profiles for the x dependence 
lead to the same spin dynamics. 

The above analytical regime no longer holds when the 
transverse motion couples with the spin through the term 
XnPxCTy in Eq. jT)). Complete analytical solutions for the 
Hamiltonian including the full spin-orbit interaction are 
not available and numerical calculations are needed to 
solve the problem. The first question arising is how the 
band structure of the simplified Hamiltonian is modified 
when transversal spin-orbit coupling is included. Figure 3 
shows the corresponding band structure for the same wire 
parameters considered above. In order to clarify the effect 
of the transversal coupling the longitudinal kinetic term 
h^ky/2m* has been subtracted and the analytical energies 
are shown with dotted lines. It can be seen that for low 



values of ky there is good agreement between the simpli- 
fied and full Hamiltonians; but this is not true when we 
approach a region of band crossings. When the {n + l)-th 
down subband reaches the up branch of the n-th subband, 
the transversal coupling produces a large effect, removing 
the band degeneracy and opening a gap between these 
two branches, (n -I- 1,|) and (n, t) 0- The gap opened 
clearly depends on the transversal subbands involved; the 
higher transversal state the larger the gap, as easily no- 
ticed from Fig. 3. It can also be observed that the range of 
ky^s where there is a good agreement between full calcu- 
lations and analytical model disminishes as the subband 
index increases. These changes on the level structure sug- 
gest significant modifications of the spin precession due to 
the coupling with transversal motion. 

As shown by Fig. 3, transversal coupling effects maxi- 
mize at the analytical model crossings, occurring approx- 
imately when 2Xf(ky — hiOQ. The fulfillment of this con- 
dition establishes a rough criterion to quantify the rele- 
vance of the transversal spin-orbit coupling. Nevertheless, 
it must be also stressed that the coupling depends sensi- 
tively on the subband index n. 

The influence of the transversal coupling can also be 
observed in the structure of the eigenspinors. In the an- 
alytical model the eigenfunctions are states having well- 
defined spin orientation. The transversal spin-orbit cou- 
pling mixes the different subbands leading to a more com- 
plex structure of the eigenstates. Figure 4 represents the 
densities corresponding to three lowest eigenstates in the 
region of maximum band anticrossing, i.e., ky = 1.70 x 10^ 
cm^^. The lowest energy state corresponds to a branch 
without crossings and, therefore, its eigenvalue and eigen- 
spinor are very similar to those of the analytical model; 
spin density mainly oriented in —x direction with a resid- 
ual structure in z orientation. The other two eigenspinors 
differ substantially from those of the analytical model; 
spin has no predominant orientation and displays a com- 
plex texture very similar to those found in Ref. ^U]. Charge 
and spin densities are related to the eigenspinor compo- 
nents through the usual expressions 



p{x) = |0„fc^-r(a:)p + \(l)nkyiix)\'^ 
m^{x) = \(f)nky'\(x)\^ - |(/'nfe„i(a;)P 

ma;(a;) = 2Re|(/)„fc„T(2;)Cfc„x(a;)} . 



(15) 



Until now we have presented the effects produced by 
the full spin-orbit coupling on the level scheme with re- 
spect to the analytical model, having results in good agree- 
ment with those of Ref. [Sj. The question that we address 
next is how this modified band structure reflects on the 
spin dynamics. We consider flrst the case of a particle 
having well-defined spin orientation at initial time and 
y-momentum in the range of band anticrossings, i.e., an 
spinor as in (|12|l evolving with ky — 1.70 x 10^ cm^^ in 
Eq. IjlOl) . The initial condition thus reads 



fky{x,0) = 



1 



ri/4 



m Wo 



1/4 



(16) 
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The spatial part of this injected spinor is equal to that 
of the lowest analytic subband. In this way, we ensure 
that it is in the low energy range, very close to the first 
transversal subband. 

The dynamical evolution of this injected particle is 
shown in Fig. 5. The spin trajectory on the Sy-Sz plane 
is quite involved, manifesting a multifrequential evolution. 
Indeed, the Fourier transform of the series shows that spin 
evolves in time as a superposition of several frequencies 
corresponding to the energy transitions between the dif- 
ferent anticrossed branches. Since the gap opened by the 
transversal spin-orbit coupling varies from one subband 
to another, several precession frequencies are found in the 
anticrossing region. On the contrary, the analytical model 
predicts a common transition energy for all the transversal 
subbands (see Fig. 3), leading to a unique precession fre- 
quency. This feature proves that spin evolution in 2D wires 
noticeably depends on the transversal state of the particle. 
In the particular case of Fig. 5 the precession is character- 
ized by two dominant frequencies, corresponding to the 
transitions between the lowest branch and the next two 
higher branches, represented by dark-gray arrows in Fig. 
3. This two-frequency evolution yields a variable spin am- 
plitude, in agreement with the displayed trajectory. The 
spin-orbit coupling makes the precessional frequencies de- 
pend on ky in a non-trivial way. At the same time, it gives 
rise to a /c^-dependent spin modulation angle that would 
lead to a certain degree of decoherence in the spin trans- 
port through the wire. 

Figure 6 shows the precessional spectra corresponding 
to an initial spinor having well-defined spin (in +z di- 
rection) and different transversal states. The precessional 
spectrum corresponding to the first transversal subband 
{n = 0) is that represented in Fig. 5. The precessional 
spectra for the higher transversal states, n = 1 and n = 
2, exhibit an enhanced multi-frequential character, since 
more transitions significantly contribute to the spin evo- 
lution. It can also be seen that the strength spreads over 
a higher frequency interval as the subband index n is in- 
creased. For low fcj,'s, out of the anticrossing region, the 
spin evolution recovers a single- frequency behaviour, with 
slight deviations from the analytical results (cf. Fig. 2). 
Note that, in general, a pure up {+z) spin injected into 
the wire will couple with the different spinor eigenstates 
and, therefore, the multifrequential evolution can not be 
avoided. 

Another question of relevance is how the transversal 
size influences the precession. From Fig. 4 we can estimate 
the width of the wire in 90 — 120 nm for hwo = 3.5 meV. 
If we reduce the value of ujq and maintain the intensity 
of spin-orbit coupling the crossings of transversal bands 
become more frequent as subbands are less spaced. Fig- 
ure 7 shows the corresponding band structure in the case 
hcoo = 0.70 meV; from the densities of the eigenspinors an 
estimated width of 200 — 250 nm can be obtained for this 
wire. Deviations from the analytical band structure are 
clearly visible in spite of the reduced ky range displayed. 
This feature is reflected in the precessional spectra which 
are fragmented even for relatively low ky^s. By contrast. 



a wire of Kluo — 14 meV, having an estimated width of 
40 — 60 nm, shows a band structure very close to that of 
the analytical decoupled Hamiltonian, as displayed in Fig. 
8. This also reflects on the precessional spectra, which are 
single-moded and match the value given by the analytical 
formula Tiujp = 2X}iky. 

In practical implementations of the Datta and Das spin 
transistor 3 the above discussed spin decoherence at sub- 
band crossings would strongly affect the device operation 
and, therefore, it should be carefully considered in the de- 
sign of appropriate working parameters. We can quantify 
the decoherence by representing the spin amplitude after 
the first spin fiip, i.e., after the first tt rotation. Figure 9 
represents this quantity as a function of wire width (loq) 
for a ky value in the anticrossing region, « huJo/2Xi^. It 
can be seen that as the wire becomes wider (lower ujq), 
the multi-frequential precession leads to an important re- 
duction of the z spin after the first inversion. 

4 Summary 

In summary, we have studied the electronic spin preces- 
sion in quantum wires with Rashba spin-orbit coupling by 
means of real time simulations. A detailed study of the 
band structure for this system has been presented and 
compared with a decoupled analytical model. It has been 
shown that deviations from the simple spin precession 
regime are determined by the relative importance of the 
longitudinal spin-orbit energy (A^fcy) and the transversal 
confining energy (huo) and that the spin precession frag- 
ments into several frequencies whose weights depend on 
the transversal state of the particle. The dependence of 
the spin precession on the wire width has also been stud- 
ied and illustrated through the band structure. 

This work was supported by Grant No. BFM2002- 
03241 from DGI (Spain). 
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Fig. 1. Fast Fourier transform (FFT) in logaxithmic scale cor- 
responding to the oscillations of a injected wavepacket evolving 
with the simplified Hamiltonian (without transversal spin-orbit 
coupling). Dotted lines indicate the analytical eigenvalues la- 
beled by the transversal subband index n and spin orentation 
s. Right panels show the non-zero densities corresponding to 
the two lowest eigenstates; charge and spin densities have been 
shifted in panel b to better distinguish them. 



Fig. 2. Lower panel: time-evolution of the spin expectation 
values in the Sy — Sz plane. Upper-left panel: Fourier spectrum 

(logarithmic scale) corresponding to the {az){t) time series. 
The dotted line gives the analytical energy. Upper-right panel: 
trajectory in the Sy — Sz plane of the {a) vector. 
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Fig. 3. Band structure, excluding the constant kinetic term 

corresponding to the first five transversal subbands 
of a wire with hoJo — 3.5 meV. Results for the full Hamilto- 
nian (solid lines) and for the decoupled one (dotted lines) arc 
shown. Transitions from the lowest down branch in the decou- 
pled (light grey) and full (dark grey) models are indicated by 
arrows. 
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Fig. 5. Lower panel: time-evolution of the (ay) and (ctz) ex- 
pectation values for an injected wave-packet with well-defined 

initial z spin orientation using the full Hamiltonian. Upper-left 
panel: Fast Fourier transform (linear scale) of the {az){t) se- 
ries; the dotted line indicates the precession frequency given 
by fiu!p = 2XRky. Upper-right panel: trajectory in the Sy-Sz 
plane of the (a) vector. 
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Fig. 4. Charge and spin transversal profiles for the three low- 
est eigenstates of the full Hamiltonian in the anticrossing re- 
gion. The panel labels correspond to the a, b and c branches 
indicated in Fig. 3. 
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Fig. 6. Precessional spectra obtained from the time-evolution 
of harmonic spinors with well-defined z spin corresponding to 
the n = 1 and 2 transversal subbands in the anticrossing re- 
gion. Dotted lines indicate the energy hu)p. 
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Same as Fig. 7 but for huio = 14 meV (Note the en- 
:;,-scale with respect to Fig. 7). 
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Fig. 7. Band structure without the longitudinal kinetic term 
h^ky/2Tri' corresponding to the first five transversal subbands 
of a wire with Tiujq — 0.7 mcV for the full Ifamiltonian (solid 
lines) and for the decoupled one (dotted lines). 
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Fig. 9. Amplitude of the {crz){t) signal after the first spin flip 
as a function of the transversal confinement strength fiuio- The 
inset illustrates the displayed quantity for the case fiujo = 3.5 
meV. 



